IPPP/12/67, DCPT/12/134 



Partially (in)visible Higgs decays at the LHC 

Christoph Englert,' Michael Spannowsky,' and Chris Wymant' 

'institute for Particle Physics Phenomenology, Department of Physics, 
Durham University, DHl 3LE, United Kingdom 

Both Atlas and CMS have reported a discovery of a Standard Model-like Higgs boson H of mass around 
125 GeV. Consistency with the Standard Model implies the non-observation of non-SM like decay modes of 
the newly discovered particle. Sensitivity to such decay modes, especially when they involve partially invisible 
final states is currently beyond scrutiny of the LHC. We systematically study such decay channels in the form 
of — !> AA jets + missing energy, with A a light scalar or scalar, and analyze to what extent these exotic 
branching fractions can be constrained by direct measurements at the LHC. While the analysis is challenging, 
constraints as good as BR < 10% can be obtained. 



I. INTRODUCTION 

Recent results obtained at the Large Hadron Collider by the 
Atlas and CMS experiments lli 0] have revealed the existence 
of a light Higgs candidate ||3[] with a mass of ~ 125 GeV. 
The observation of this new particle combines evidence in the 
H ^ -fj, H ^ ZZ and H W+W^ channels. Given 
the yet small collected luminosity, the properties of this newly 
discovered resonance are still subject to large statistical uncer- 
tainties 101 ■ Analyses targeting e.g. spin and CV IH^l of the 
new particle and a more precise extraction of its couplings to 
known matter will be addressed with a larger data sample. 

The current observations leave open a plurality 
of phenomenological possibilities of Higgs sector- 
modifications 10, Si- Especially the extraction of the 
resonance's couplings is influenced by non-standard de- 
cays jgl- fisll since it is based on a fit to combinations of 
various production p and decay modes d. These are functions 
of the partial and total decay widths and all couplings {gi}: 
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The total width Ftot in current fits is typically approximated 
by including a freely flowing invisible partial width fl4ll to 
the list of decay modes or by imposing the constraint Finvis ~ 

5invi,s = 0. 

Extracting such an invisible* or partially visible branching 
ratio is an experimentally ambitious task. The decay of the 
Higgs via a light scalar or pseudo-scalar A can be buried in 
a large hadronic background, experimental systematics can 
limit the sensitivity to such decays, and non-SM phenomenol- 
ogy can easily be missed. The signature H — > AA occurs in 
many extensions of the SM, e^., the next-to-minimal super- 
symmetric Standard Model lisll . Higgs-portal models lUol . 
and whenever an approximate symmetry of the Higgs poten- 
tial is explicitly broken by a small term in the potential, giving 
a light pseudo-Nambu-Goldstone boson (see IitIi for a clear 



*Note that from Eq. (T) "invisible" also means fully visible in a non-standard 
search channel. 



discussion). However the signature can be missed by standard 
searches, depending on how A itself decays, and new dedi- 
cated strategies need to be devised. Obviously, an observation 
of such a novel decay channel H — > AA would directly imply 
physics beyond the Standard Model. 

Subjet methods (pioneered in Ref. ifisll ) have proven par- 
ticularly successful in getting a handle on such a modified 
phenomenology. In particular, subjet analyses applied to the 
decay chains H — > 2 A — > AX have unravelled potential sen- 
sitivity to these non-standard decays if the (pseudo)scalar A 
is hght 0(10 - 20 GeV). A decay of the Higgs to AA is 
well motivated on general grounds; one can keep a reasonably 
open mind with regards to how the AA subsequently decay. 
H AA Ag is considered in fH], H A A 4r 
in il, H 2A^ 2t2h in Oil and i? -» 2A 4c in i2lll. 
Higgs decays into resonances with masses close to hadronic 
bound states have been studied in Refs. ||TC)1|22| |. Sensitivity 
to the signatures discussed in Refs. Il0lll7l]l9 - l22ll follows 
from fundamentally distinct QCD and electroweak properties, 
highlighting the diverse power that jet substructure-based ap- 
proaches offer 

There is one major difference in the analysis of purely 
hadronic final states compared to electroweak final states. 
Quite often the latter involve a significant amount of miss- 
ing energy, which is aligned with the direction of the fat jet 
that is input to grooming and/or tagging algorithms (23] . Nor- 
mally such an event topology is avoided to minimize system- 
atic uncertainties. In, e.g., searches for supersymmetry in the 
jetsH-missing energy channel ll24ll one requires a missing en- 
ergy vector well isolated from a number of hard jets, to 
reduce systematics. Decays H AA with niH/mA S> 2, 
on the other hand, naturally involve non-isolation of 
which might even be not too large depending on the decay 
of A. In the phase space region where we can separate signal 
from background a SUSY-inspired search strategy based on a 
strong isolation becomes impossible. 

Strategies to deal with non-SM Higgs decays involving 
both missing and hadronic energy 112511 are, hence, limited. 
A direct generalization of the standard invisible Higgs decay 
search in the weak boson fusion channel iflTll is not possible. 
Such an analysis relies on a central jet veto to obtain a suf- 
ficiently large signal-over-background ratio S/B to study the 
jets' azimuthal angle correlation. Applying a central jet veto 
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FIG. 1: Higgs decay topologies in tiie simplified models that we study for the purpose of this paper, "vis" is a placeholder for uu, dd flavor 
quark pairs that give rise to visible hadronic energy. 



to central decay H — > ^-p+hadrons removes the signal that 
we would like to investigate. 

Adapted Higgs+monojet searches {^[iltl are challenged by 
overwhelmingly large dijet and weak boson+jet backgrounds 
and trigger issues, as soon as missing energy from invisible 
Higgs decays decreases when turning to a partially visible de- 
cay. The latter can at least be partially cured by focussing 
on the Higgs' pT-distribution's tail. This, however, comes at 
the price of larger theoretical and experimental uncertainties, 
which imperatively need to be included to reach a realistic 
formulation of the branching ratio constraints. 

This leaves Higgs-strahlung pp HZ as the best- 
motivated channel to study the situation we have in mind at 
the LHC. Trigger issues are avoided by reconstructing the lep- 
tonic Z e^e^ , fi^i^i^ decay and no further adjustments to 
the trigger settings or thresholds are necessary to perform the 
measurement(this process is also studied in the context of in- 
visible Higgs searches ll26ll ). Furthermore, the jet energy scale 
calibration is performed with Z+jet events, whose distribu- 
tion is both theoretically and experimentally under good con- 
trol ll27ll28ll . Data-driven methods can be straightforwardly 
applied in "ABCD" approaches, when e.g. comparing boosted 
to the un-boosted Z boson kinematics. 



II. DECAY TOPOLOGIES 

We employ a simplified-model based approach ll29ll to in- 
vestigate the LHC search potential to the H ^y+hadrons 
signature. More precisely we study the Higgs decay realiza- 
tions depicted in Fig. [T] These scenarios (a)-(d) are charac- 
terized by different kinematics and different relative contribu- 
tions of missing and hadronic energy in the final state. 

We limit our analysis to the light flavor final states u, d; 
particles with ~ 10 GeV masses and a significant coupling 
to b quarks are not only constrained by upsilon measure- 
ments llsoll . but also give rise to plethora of dedicated phe- 
nomenological handles on the final state 1221 [31II32I1 that we 
do not wish to exploit to be as general as possible. Similar 
arguments hold for the H ^ At decay, which gives rise to 
sparse but focused hadronic energy deposits ||20|]. Both these 
cases (and also H — > 4g) are fundamentally different from 
the topologies of Fig. [T] where "visible" refers to light flavor 
quarks. The quarks undergo normal showering and hadroniza- 
tion leaving neither the possibility for flavor tags nor for de- 
tecting pronged decays from counting charged tracks, as done 
in the h- and r-flavored decays of H ^ AA, respectively. 
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FIG. 2: A$rp/frp as a function of the true $rp that follows from 
particle flow and early LHC data (ssll . We also display the fitted 
function that is employed for our analysis. 



III. ELEMENTS OF THE ANALYSIS 

We implement the decay topologies using Feyn- 
RULES lUsj] and use its interface to Sherpa to gen- 

erate events for the mass choices mn = 125 GeV, niA = 
20 GeV, TOinvis = 10 GeV. This choice is not special and 
the details of eventually extracting the branching ratio is not 
sensitive to the particular value unless niH iriA- We 
will comment on the possibility to extract niA in Sec. lIVI 

We generate background events using Sherpa and include 
M^VF-njet, WZ+jet, ZZ+jet, and tt+jets as the main back- 
grounds to om pp -+ {Z ^ t^l^) + jet + firj. analysis. 
We normalize our signal and background event samples to the 
corresponding higher order-corrected cross sections Isil - fssll . 
Studying the impact of a mismeasurement of Z-njets events re- 
quires the simulation of a realistic detector environment, and 
should be addressed by the experiments. However, we may 
assume on the basis of Refs. ll27il28ll that this background can 
be brought under sufficient control and can be subtracted from 
the eventual distribution also when the missing energy vector 
is collinear to the jet. We include a flat shape uncertainty of 
the background distribution to partly account for the jet en- 
ergy scale uncertainty in the computation of the expected BR 
Umits in Sec.HVl 

Associated Higgs production with SM-like Higgs decays to 
66 and rr also comprise backgrounds to our (Z — > t^l^ ) + 
^j-+jet analysis and we include them consistently through- 
out. 
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FIG. 3; Missing transverse energy (a) and fat jet transverse momentum (b) of the decay scenarios Fig.[TJa)-(d) and the contributing backgrounds 
after all analysis steps have been applied. 



We reconstruct the events' visible final states using a hybrid 
ECAL+HCAL implementation which granularizes the final 
state particles on grids with Arj x Acj) = 0.025 x 0.025 (0.1 x 
0.1) ECAL (HCAL) as massless cell entries and feed the re- 
constructed objects to a smearing routine which mimics detec- 
tor effects as described later on. Doing so, we reconstruct the 
full three-momenta from the detector geometry with an invari- 
ant mass = for each ECALh-HCAL cell, which contains 
a single or multiple hits. 

In the analysis, we first reconstruct isolated stable lep- 
tons, by requiring the hadronic energy deposit in the vicin- 
ity of the lepton candidate (AR = ^/Aj]"^ + Acj)^ < 0.1, 
where 77, and cj) are pseudorapidity and azimuthal angle respec- 
tively) to be smaller than 10% of the lepton candidate's trans- 
verse momentum. We furthermore require exactly two leptons 
(PT.e > 10 GeV, \T]e\ < 2.5) of identical flavor and oppo- 
site charge that recombine the Z mass within mz ± 10 GeV 
(note that this also removes higher-pronged H tt events) 
and demand pt,z = [pii + P12)t > 130 GeV. Subse- 
quently, we cluster Cambridge/Aachen jets with i? = 1.5 us- 
ing FastJet ||40Ii and we require at least one such fat jet with 
Pt > lOO GeV. We apply the mass drop and asymmetry 
criteria to this jet as explained in Ref. ifisll . we filter and trim 
it ll4lll . and we keep the three hardest subjets ifisll . Further- 
more, the jet needs to pass a 6-veto in |r/| < 2.5, for which we 
choose a working point ll42ll with a flat 70% tagging efficiency 
and a flat 10% fake rate. The b veto removes the pollution 
from H, Z ^ bb and significantly reduces the tt background. 

In the next step we remove the fat jet from the event and 
re-cluster the remaining tracks to anti-fc^ jets with R = 0.4 
and Pt > 30 GeV. In order to project out the signal region 
we are interested in, we require a large missing transverse en- 
ergy Et > 50 GeV, which is reconstructed from all visible 
hits. It needs to be isolated from the reconstructed Z boson 
R{$Tj Z) > 2, and we veto additional jets around the Z: 
Vj : R{Z, j) > 1.5 to further reduce events. 

Throughout the analysis we include detector smearing 



which we choose for jets and leptons according to Ref. 11431] : 
AE 5.2 0.16 



jets 



AE 

leptons : —— = 0.02 , 

E 



(2) 



and we include the missing energy response from recent par- 
ticle flow fits of CMS iH (obtained from flie 7 TeV data set, 
Fig.|2j: 



missing energy 



Ati 



2.92 



0.07. 



(3) 



In fact, the quoted jet and lepton energy resolutions are better 
than the ones obtained from particle flow observables ll39tl . 
The resolution of supersedes the original estimates of 
Ref. ll43ll and we can therefore assume that Eq. (O is a con- 
servative baseline for the 14 TeV run at luminosities at which 
the limits on the BRs will be stringent. We note that the fx 
resolution is the main limiting factor of the analysis that we 
pursue. 

After applying all cuts, we are left with only minor possibil- 
ities to enhance signal over background (apart from including 
an invariant mass cut around the reconstructed Higgs mass, 
of course). Fig. [3] displays the $Jrp and fat jet pT distribution 
for the signal (assuming BR = 1 to the discussed topologies) 
and the discussed backgrounds. With all described cuts ap- 
plied, tt is sufficiently reduced, while the electroweak yy-njet 
backgrounds are still large. Additionally, the fj- ™d px dis- 
tributions do not show any particular discriminative features, 
which would allow to gain in S/B from tighter and harder 
missing energy and jet pT cuts (note that this especially holds 
when we need to relax BR = 1 for deriving upper limits on 
the branching ratio). We find that SM-like Higgs decays are 
sufficiently suppressed to make negligible contribution to the 
background. 

In order to derive expected limits on the branching ratios 
we need to reconstruct the 125 GeV resonance in this search 
channel in the most efficient way. There are, in fact, dedicated 
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FIG. 4: Transverse mass or transverse cluster mass distributions after filtering and filtering+trimming, depending on the scenario of Fig.[T] We 
assume BR(scenario (i)) = 1, i = (a) — (d) for illustration purposes. 



observables, which faciUtate the reconstruction of invariant 
mass peaks when missing energy is involved in the decay of 
a heavy resonance and they also prove helpful in the present 
context. One example which is useful for our purposes is the 
transverse mass definition 




which generalizes the invariant mass definition of Ref. ll44ll to 
the present situation. In Ref. jiill an analogue definition of 
the invariant mass is applied to i7 —> tt (see also lf32ll ) with 
full leptonic t decays in the collinear approximation, i.e. the 
individual |>-lepton correlation from each r decay is incorpo- 
rated. Applied to our signatures this corresponds to a leptonic 
version of Fig.[Tla), and we can expect a good reconstruction 
of the Higgs resonance for this decay scenario. 
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FIG. 5: 95% confidence level exclusion of the various scenario's branching ratio (or signal strength ^ — a x BR/itsm to be more specific) 
from a direct measurement along the lines of Sec. [in] We assume a background template uncertainty of initially 10%, which is we assume to 
saturated at 5% at end-of-lifetime LHC luminosities of £. = 1000/fb. 



Another observable which reconstructs the mass of a heavy 
decaying resonance from a kinematic endpoint measurement 
is the so-called transverse cluster mass 14511 

mT,c = +^7^) -{pT..3+pTy- (5) 

This observable is a good choice if there is a relatively 
large amount of visible energy compared to invisible en- 
ergy ll20i,|4^ . as encountered in our scenarios Fig. [TJc) and 
(d), for which it allows to reconstruct the Higgs mass from a 
Jacobian peak. 

For signatures involving a lot of missing energy such as 
Fig. [He), we need to break the degeneracy of the $x vector 
to facilitate an approximate mass reconstruction. For scenario 
Fig. nib) we project p onto the hardest reconstructed subjet 
after filtering j, and assign an energy $ = E-^ We define 
a second missing energy four vector from the shift with re- 
spect to the original missing energy vector ^ = |>^^.^ — p. This 
breaks the missing energy degeneracy and assumes a collinear 
democratic A decay in the lab frame. Resolution effects have 
the biggest influence on this decay scenario and the mass re- 
construction will be worst among the four scenarios that we 
discuss in this work. 



The corresponding invariant mass distributions for our sce- 
narios Fig.[na)-(d) after all cuts have been applied are shown 
in Fig.m where we again assume BR(scenario (i)) = 1, i = 
(a) — [d). We pick the best representation of the Higgs 
mass peak by studying the invariant mass distributions be- 
fore and after filtering and trimming. From the different jet- 
substructure that follows as a consequence of the different de- 
cay topologies, trimming can remove too much of the signal 
when the light final state quark pairs are relatively widely sep- 
arated in scenario (a). 



When focussing on scenarios with an increasingly large 
amount of $j, compared to hadronic activity, sensitivity de- 
creases. This is a consequence of the applied subjet analysis, 
which requires by definition a significant amount of a hadronic 
energy. Depending on the scenario we therefore become in- 
creasingly sensitive to uncorrelated initial state radiation and 
underlying event splash-in which degrades the mass resolu- 
tion, while filtering and trimming serve to suppress part of the 
pollution from underlying event. 



6 



0.6 




filtered+trimmed mx,c [GeV] 
(a) 



FIG. 6: (a) Transverse cluster mass for scenario (c) with the additional requirement Eq. Q, and (b) the resulting branching ratio extraction 
along the lines of Sec. lIVI 



IV. BRANCHING RATIO LIMITS AT THE LHC 

We move forward and use the distributions of Fig.|4]to de- 
rive an estimate of the bounds on the branching ratios that can 
be obtained at the LHC 14 Te V from direct measurements. We 
follow the experimentalists' approach and use the histograms 
of Fig. m as an input to the CL^ method HS] to derive 
95% CLs exclusion limits on the branching ratio using the 
binned log-likelihood test statistic. We include a shape uncer- 
tainty of the background template in the computation of con- 
fidence levels that eventually enter the CL^ ratio. We choose 
this shape uncertainty as a flat profile which we assume de- 
pendent on the integrated luminosity 

A(£) = 10% - 5%y^C/1000fb-^ (6) 

to mimic an improved understanding of the measurement and 
theoretical uncertainty in due time. 

The expected 95% confidence level exclusion is shown in 
Fig. |5] Direct measurements yield results BR < 0.1, . . . , 1, 
and are highly sensitive to the amount of hadronic energy that 
we observe in the final state and the correlation of fx with 
the visible part of the final state. 

We learn from Fig. |4(b)| that the combination of little 
hadronic energy, trigger criteria and detector effects together 
with initial-state-radiation pollution results in an extremely 
challenging signature at the hadronically busy LHC environ- 
ment. As a direct consequence, only weak limits can be ob- 
tained for decay scenario (b). This is mostly due to the very 
limited possibilities to improve S/ B and the poor mass res- 



olution that can be obtained from a signature that predomi- 
nantly governed by missing energy in a channel whose dibo- 
son backgrounds' also contain a sizable amount of 

Turning to final states that contain more hadronic energy 
Fig. [TJc) and (d), we can impose stronger limits from direct 
measurements (BR < 0.3). An optimal choice of observables 
adopted to the specifics of Fig. [Tta) in comparison to (c) can 
push the bounds to the 10% level. 

An alternative route for putting limits on scenario (c), where 
hadronic energy correlations in the fat jet's substructure can 
be resolved, is exploiting the fat jet's "active area" (see also 
Ref. ll20ll ) in a way that also incorporates the correlation with 
the missing energy. This is straightforwardly achieved by im- 
posing an additional cut on the ratio of the fat jet's trans- 
verse momentum and the events reconstructed transverse clus- 
ter mass mT.c{$, j) 

PTj/mT,c>2. (7) 

The resulting rriT.c is shown in Fig. |6] together with the pro- 
jected BR extraction. While the reduction in statistics eventu- 
ally compensates the enhancement in S/ B (compare Fig.|6jb) 
to Fig. He)), we are able to approximately reconstruct our 
light Higgs partner = 20 GeV from the jet substructure. 
This is only feasible by exploiting visible final state energy 
coiTelations - identical strategies proof unsuitable for scenar- 
ios (a), (b) and (d). 

A successful and intangible extraction of the light scalar 
mass from the associated production channel hence depends 
on a good understanding of the background at small tot.c and 
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a sizable branching ratio H — > AA — > scenario (c). Current 
fits limit the invisible branching ratio to be BR < 0.5. Since 
all of our signatures are missed in standard searches, our sce- 
narios are constrained by this loose bound, so that there is the 
possibility to find light scalars in that particular mass range in 
the near future. 



V. SUMMARY AND CONCLUSIONS 

After the discovery of the Higgs-like resonance at around 
125 GeV, a further investigation of the resonance's com- 
patibility with the SM expectation is ongoing. A precise 
measurement of the resonance's couplings to SM matter de- 
pends on the treatment of non-standard partial decay widths, 
which is typically performed in global fit to the quantita- 
tive (non-)observation in different exclusive Higgs final states. 
Non-standard situations are hereby often treated as a nuisance. 

Hence, the necessity to investigate non-SM phenomenolog- 
ical situations is imperative: If an indication of a non-zero 
branching ratio in a specific channel is accessed in a direct 
measurement, the evidence for physics beyond the SM will be 
accompanied by an insight in its particular realization as op- 
posed to a mere parametrization of Finv in Eq. ^ as part of 
a global fit. This inevitably comes at the price of a decreased 
sensitivity (i.e. less stringent limits) which results from in- 
vestigating a single channel as opposed to combinations and 
correlations of many. 

We have investigated the possibility to constrain non- 
standard branching ratios via direct measurements for a broad 
range of effective models, which contain both invisible and 



visible final state particles in this paper. Gaining sensitiv- 
ity to these phenomenological situations is characterized by a 
number of shortcomings, that range from trigger issues all the 
way to the Higgs mass reconstruction at a reasonable signal- 
over-background ratio. Since trigger thresholds are typically 
bound to hardware specifications, we devise an analysis strat- 
egy which builds upon associated production and a boosted 
Higgs final state. As a consequence, we need to deal with 
an event topology where missing energy is collimated with 
the visible Higgs decay products in decays H — > AA with 
2to^ ^ 125 GeV (a choice 2m^ ~ 125 GeV exhibits too 
much ^j^-jet decorrelation to facilitate signal reconstruction 
on top of a large background rejection). Recent investigations 
targeting the reconstruction of based on particle-flow as 
well as the extraction of fake from a combination of the- 
ory and experiment suggest that systematics can be brought 
under sufficient control to perform such analyses in this par- 
ticular channel. 

We find that, depending on the specific realization of the 
partially visible Higgs decay (which is characterized by dif- 
ferent ^y-jet substructure correlations) bounds as good as 
BR < 10% can be obtained from a measurement of boosted 
Higgs final states, applying a rather generic cut set up which 
serves to obtain a reasonable S/B. Constraints can subse- 
quently be obtained from an optimal choice of observable 
which, depending on the amount of visible final state energy, 
also allows to partially reconstruct parts of the Higgs' poten- 
tial decay chain. 
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